Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 12 September 2012 (MN I^TeX style file v2.2) 



A suppressed contribution of low mass galaxies to 
reionization due to supernova feedback 

J. Stuart B. Wyithe 1 , Abraham Loeb 2 

1 School of Physics, University of Melbourne, Parkville, Victoria, Australia 

2 Astronomy Department, Harvard University, 60 Garden Street, Cambridge, MA 02138, USA 
Email: swyithe@unimelb.edu.au; aloeb@cfa.harvard.edu 



12 September 2012 



ABSTRACT 

Motivated by recent observations of the star formation rate density function out to 
z ~ 7, we describe a simple model for the star formation rate density function at 
high redshift based on the extended Press-Schcchter formalism. This model postulates 
a starburst following each major merger, lasting for a time igp and converting at 
most /* max of galactic gas into stars. We include a simple physical prescription for 
supernovae feedback that suppresses star formation in low mass galaxies. Constraining 
£sf and /*. ma x to describe the observed star formation rate density at high redshifts, 
we find that individual starbursts were terminated after a time tsF ~ 10 7 years. This 
is comparable to the main-sequence lifetimes of supernova progenitors, indicating that 
high redshift starbursts are quenched once supernovae feedback had time to develop. 
High redshift galaxies convert ~ 10% of their mass into stars for galaxies with star 
formation rates above ~ 1M© yr _1 , but a smaller fraction for lower luminosity galaxies. 
Our best fit model successfully predicts the observed relation between star formation 
rate and stellar mass at z >4, while our deduced relation between stellar mass and 
halo mass is also consistent with data on the dwarf satellites of the Milky Way. We 
find that supernovae feedback lowers the efficiency of star formation in the lowest 
mass galaxies and makes their contribution to reionization small. As a result, photo- 
ionization feedback on low mass galaxy formation does not significantly affect the 
reionization history. Using a semi-analytic model for the reionization history, we infer 
that approximately half of the ionizing photons needed to complete reionization have 
already been observed in star-forming galaxies. 
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1 INTRODUCTION 

The galaxy luminosity function is the primary observable 
that must be reproduced by any successful model of galaxy 
formation. At z J>6, it also represents one of the most im- 
portant observables for studying the reionization of cosmic 
hydrogen. Developing a theoretical picture of the important 
processes involved in setting the star formation rate at high 
redshift lies at the forefront of understanding this important 

|Finlator et al.|[20TT| 



et al. 20111. Observations show that the value of de- 



2010 



cosmic epoch (e.g. Trenti et al. 
Mu&oz fc Loeb||2011||Raicevic et al.||2011| |Salvaterra et al.| 
20111). 



The luminosity function of Lyman-break galaxy candi- 
dates discovered at z > 6 in the Hubble Ultra-Deep Field is 
described by a Schechter function with characteristic den- 
sity 'J* in comoving Mpc -3 , and a power-law slope a at lu- 



minosities L below a characteristic break L* (e.g. Bouwens 



creases towards higher redshift as expected from the dark 
matter halo mass-function (e.g. Mufioz fc Loeb|2011 1, while 
the faint end slope of a « —1.8 is observed to be roughly 
independent of redshift ( McLure et al.|20'09 Bouwens et al. 
|2011[ ). A complication that arises when modelling the lumi- 
nosity function is that models predict a star formation rate, 
which must then be converted to a luminosity assuming an 
initial mass function (IMF) for the stars. While this calcu- 
lation is straightforward, of more importance is the poten- 
tial contribution of reddening. As shown by Bouwens ct al. 
(20121, the blue continuum slope of z J>6 galaxies is de- 
pendent on luminosity and redshift. As a result there is a 
dust correction which is both luminosity and redshift de- 
pendent and difficult to reproduce from first principles re- 
liably in a simple model. This renders theoretical studies 
of the luminosity function difficult to interpret. For exam- 
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pie, in the model of Mufioz & Loeb (2011 1 a distribution of 



mass-to-light ratios is assumed whose mean is independent 



of mass, while in Trenti et al. (20101 the mass-to-light ratio 



is assumed to be independent of redshift. Neither of these 
assumptions holds based on the more recent observational 
work of Bouwens et al. (20121. More complex models are 



able to better reproduce many of the observed properties, 
and so make more robust predictions of the physics. For ex- 



ample Finlator et al. (20111 have modelled the growth of 
stellar mass in high redshift galaxies using hydrodynamical 
simulations coupled with sub-grid models for processes in- 
cluding star formation and metal enrichment, and broadly 
reproduce the luminosity function evolution as well as the 
blue colours of the young stellar populations at high redshift. 
Similarly, |Salvaterra et al.| ( |2011[ ) and |Jaacks et al.| ( |2012[ ) 
have calculated the evolution of the luminosity function in 
detailed numerical simulations including calculations of en- 
richment and dust reddening, with the latter also including 
additional physics related to the transition from population- 
Ill to population-II stars. A focus of these numerical studies 
is the role of supernova (SNe) feedback on star formation, 
particularly in low mass systems. In a separate approach, 



Raicevic et al. (20111 (see also Benson et al. 2006 Lacey 



et al.||2011[ ) used a semi-analytical galaxy formation code 



based on Monte-Carlo merger trees to model the evolution 
of the high redshift luminosity function, and study the ef- 
fect of SNe feedback on the global ionizing photon budget 



and global ionization. In particular jRaicevic et al. (20111 



evaluated the ionizing photon budget, finding that although 
galaxies should produce sufficient ionizing photons to com- 
plete reionization, most of the galaxies responsible would be 
below the detection threshold of current surveys. 

A simpler way to constrain theory by observations is 
to estimate the star formation rate density observationally, 
where the correction is made from luminosity to star for- 
mation rate using the observed continuum properties of the 
galaxies under study. A powerful probe of the physics of star 
formation is then provided by the star formation rate den- 
sity function (i.e. the number of galaxies per unit volume 
per unit star formation rate). Recently, this has become a 
viable approach following the work of Smit et al. ( 2012 1 who 



combined estimates of dust extinction at z ~ 4 — 7 with mea- 
surements of the UV luminosity function in order to derive 
star formation rate density (SFRD) functions at z ~ 4, 5, 
6 and 7. These SFRD functions provide a physical descrip- 
tion of the build-up of stellar mass in galaxies at high red- 
shift. The resulting SFRD functions are well-described by a 
Schechter function, with a characteristic break separating a 
shallow dependence of SFRD on star formation rate at low 
luminosities from the exponential dependence at high lumi- 
nosities. As mentioned, the physics of star formation in high 
redshift galaxies has important implications for the process 
of reionization. Of particular interest is the star formation 
rate in low mass systems. For example, it has been shown 
that the growth of HII regions during reionization may be in- 
fluenced by radiative feedback in the form of suppression of 
galaxy formation below the cosmological Jeans mass within 
a heated intergalactic medium (IGM) ( jDijkstra et al.|2004 |, 
although the importance of this effect remains controver- 



sial (Mcsinger & Dijkstra 2008). If present, this radiative 



suppression of low mass galaxy-formation delays and ex- 
tends the reionization process, which though started by low 



mass galaxies, must then be completed by relatively mas- 
sive galaxies (e.g. Iliev et al.|2007 1 . On the other hand, the 
semi-analytic modelling of Raicevic et al. (20111 indicates 



that photoionization suppression of star formation in these 
low mass galaxies is unlikely to significantly affect reion- 
ization, because of feedback effects on the star formation 
efficiency implies that their contribution to reionization is 
small. The contribution of galaxies to the reionization of 
the IGM is dependent on the star formation rate and ini- 
tial stellar mass function, but is also limited by the frac- 
tion of ionizing photons that escape their host galaxies. If 
the escape fraction is small, then star formation had to be 
very efficient at high redshift in order to reionize the Uni- 
verse. The escape fraction is therefore a critical parameter in 
studies of the connection between high redshift galaxy for- 
mation and reionization. Attempts to determine the escape 
fraction have been dominated by direct observations of rel- 
atively low redshift galaxies, and by numerical simulation. 



Both observational (e.g. |Steidel et al.||2001 


Fernandez-Soto 


et al.|2003 


Shapley et al.|2006||Siana et al. 


20071, and theo- 


retical (e.g. 


Razoumov & Sommer-Larsen|2006||Yajima et al.| 


2009 Gnedin et al.|2007 Wise & Cen|2009 1 estimates of the 



escape fraction at the Lyman-limit are currently uncertain, 
with an expected range of 0.01 <C/ CS c <T- 

In this paper we aim to utilise this new determination of 
the build up of stellar mass at high redshift to constrain star 
formation scenarios in high redshift galaxies, with particular 
attention to the possible consequences of SNe feedback for 
the reionization of hydrogen. In § [2j we describe an analytic 
model for the star formation rate density function at high 
redshift based on the extended Press-Schechter formalism 
( Lacey fc Cole|1993[ ) with a simple physical prescription for 



SNe feedback that suppresses star formation in low mass 
galaxies. We confront our model with the observed SFRD 
function in § [3] in order to constrain the starburst lifetime 
and the stellar mass - star formation rate relation, and show 
the predictions of the model for the escape of ionizing pho- 
tons from star-forming galaxies in § [4] We investigate the 
implications for the reionization history in § [5]|6j and con- 
clude in § [7] In our numerical examples, we adopt the stan- 
dard set of cosmological parameters ( Komatsu et al.|2011 |, 



with values of fi b = 0.04, Q m = 0.24 and f2 A = 0.76 for the 
density parameters of matter, baryon, and dark energy, re- 
spectively, h = 0.73, for the dimensionless Hubble constant, 
and as = 0.82. 



2 MODEL 

The star formation rate in a galaxy halo of mass M that 
turns a fraction /* of its disk mass m<jM into stars over a 
time tsF is 



SFR = O.lSMeyr" 1 



0.17 



0.1 J V1O 8 M 



tSF 

10 7 yr 



(1) 

We assume that major mergers trigger bursts of star forma- 
tion, and constrain the starburst lifetime required to repro- 
duce the observed star formation rate density function. The 
star formation rate density function (i.e. galaxies per Mpc -3 
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per unit of SFR) can be estimated as 



$(SFR) 



Eduty 



/ dN 2 
I AM t„ 



dtdAM 



Mi, AM 



dn 
dM 



f dSFR Y 
\ dM J 



(2) 



where eduty is the fraction of the Hubble time (£h) over which 
each burst lasts, and dn/dM is mass function of dark mat- 
ter halos ( |Press fc Schechter||1974| |Sheth fc Tormen||1999 1. 
The rate of major mergers (d-/V m erge/cK) is calculated as the 
number of halos per logarithm of mass AM per unit time 
tha t merg^jwith a halo of mass Mi to form a halo of mass 
M ( Lacey fc Cole|1993 |. We assign a 2:1 mass ratio to ma- 



jor mergers (i.e. Mj = |M and AM = M/3). The ionizing 
photon rate from a starburst drops rapidly once the most 
massive stars fade away on a timescale of t B ~ 3 x 10 6 years 
(Barkana & Locb 200Tj) . If the starburst lifetime isF is much 
longer than t s , then the duty cycle associated with a star- 
burst is set by isF • However, if the starburst is shorter than 
t B , the star formation remains visible for a minimum of t E . 
Thus the duty-cycle is 



Eduty 



ts + tsF 



For comparison with observations we define 
<P(SFR) = In 10 X SFR x $, 



(3) 



(4) 



which has units of Mpc per dex. 

We expect that SNe feedback will alter the fraction of 
gas in a galaxy that is turned into stars (e.g. |Dekel fc Woo| 
2003). To determine the mass and redshift dependence of /* 



in the presence of SNe we suppose that stars form with an 
efficiency /* out of the gas that collapses and cools within a 
dark matter halo and that a fraction Fsn of each supernova 
energy output, -Esn, heats the galactic gas mechanically (al- 
lowing for some losses due to cooling). The mechanical feed- 
back will halt the star formation once the cumulative energy 
returned to the gas by supernovae equals the total thermal 



energy of gas at the virial velocity of the halo (e.g. Wyithe 
& Loeb 2003b). Hence, the limiting stellar mass is set by the 
condition 

M 1 

— E SN F SN f t f d = E b = -m d Mv 2 iv . (5) 
WSN 2 

In this relation Eh is the binding energy in the halo, wsn 
is the mass in stars per supernova explosion, and the total 



stellar mass is M* = m d M /*,tot where /* 



N n 



is the total fraction of the gas that is converted to stars 
during major mergers, and -/V mcrgc is the number of major 
mergers per Hubble time. The parameters ft and /a denote 
the fraction of the SNe energy that contributes because of 
the finite timescale of the SNe feedback or the disk scale 
height being smaller than the SNe bubble. These terms are 
described in more detail below. 

The ratio between the total mass in stars and dark mat- 
ter is observed to increase with halo mass as (M*/M) oc 



1 In addition to the excursion set approach of l |Lacey fc Cole| 
|l993| ), we have also computed merger rates using the fitting for- 
mulae based on numerical simulations of |Fakhouri et aTT] ( |2010] ) . 
We find consistent results using either approach. 



M 5 for Mj, < 3 x 10 10 M Q , but is constant for larger stel- 



lar masses (Kauffmann et al. 20031. Thus, the star forma- 



tion efficiency within dwarf galaxies decreases towards low 
masses. For comparison with equation |5|, a Scalo (19981 
mass function of stars has il>sn ~ 126 M© per supernova 
and Esn = 10 51 ergs, and so we find that M* = 3 X 10 10 
Mq and v c ~ 175 km/s (the typical value observed locally; 
see e.g. |Bell fc de Jong|2001 1 implies /*,tot ~ 0.1 for a value 
of -Fsn ~ 0.5. Smaller galaxies have smaller values of /*. 
Equation (JsJ indicates that 

/* = 



0.008 

-^V mcrec 



M 
10 10 M Q 



1 + z 
10 



(/t/dF Sr 



(6) 



We utilise equation (|6| with equation Q as a function of 
the parameters tsF and /*, ma x. 

2.1 Disk structure 

The effect of SNe feedback is dependent on the conditions of 
the interstellar medium (ISM) gas. We assume that the cold 
gas (out of which stars form) occupies a self-gravitating ex- 
ponential disk, with surface mass density £(r) = £oe~ r/ ' i?d , 
where 

m d M 



So = 



2-kRY 



and R d is the scale radius 
-Rd = 



V2 



(7) 



(8) 



where rrid is the mass fraction of the disk relative to the halo 
and A ~ 0.05 is the spin parameter of the halo (Mo et al. 
|1998[ ). The virial radius of a halo with mass Mhaio is given 
by the expression 



i? vir = 0.784/i" 1 kpc 



M h 



[c(*)r 



l + z 



where < = [(fi m /^)(A c /187r 2 )], f!^s[l+ (n A /n m )(l + 
z)~ 3 ]-\ A c = 18 tt 2 + 82rf - 39d 2 and d = - 1 (see 
equations 22-25 in |Barkana fc Loeb|200l[ for more details). 
The scale height of the disk at radius r is 

7rGL(r) 

We adopt the density in the mid plane at the scale radius, 
within which half the gas is contained, as representative of 
the density of the ISM. At r = R d the scale height is 



H = 



Gm d M 
0.034 kpc 



x 



A 
O05 

'1+2 



0.17 

12 



/ M 



10 



V1O 8 M 

2 



-1/3 



10km/s 



The number density of particles in the mid plane at r 
is 

_ G(m d M) 2 
np ~ 8irm p c 2 R4 x 2.71 2 ' 

where m p is the mass density of particles. 



(10) 
= Rd 

(11) 
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We note that if the gas disk becomes stable to frag- 
mentation at a radius beyond which there is a significant 
fraction of gas by mass, then the half mass radius of the 
stellar disk may not equal the scale radius of the gas-disk. 



However, Wyithe & Loeb (2011 1 find that the disk becomes 



stable (based on Toomre-Q criterion) only at 3-4 scale radii. 



2.2 Supernova evacuation of the ISM 



Clarke & Oey 1 2002| ) presented a simple analytic model for 
the effect of supernovae on the interstellar medium which we 
apply to high redshift galaxies. In this model, clusters of N e 
SNe produce super-bubbles in the ISM with a radius R e at 
which the super-bubble comes into pressure balance with the 
ISM. This radius can be found by approximating R a as the 
radius within which the thermal energy of the ISM equals 
the mechanical energy of the SNe cluster. The timescale as- 
sociated with the evacuation of a super bubble in the ISM 
by a SNe cluster is t c = 4 x 10 7 years, corresponding to the 
lifetime of the lowest mass SNe progenitor. The evacuation 
radius for a cluster of N e SNe, each with energy output _Bsn 
within an ISM of sound speed c s is therefore 



R c — 



f 2N e E S N 
\2%m p n p cl 



(12) 



yielding 



Re = 



0.08 kpc 



10 

M 



V10 8 M Q 



10 51 erg 
l + z 



A 
O05 



0.17 



10 



(13) 



The derived value of R e ignores radiative losses of the 
super-bubble before it comes into pressure equilibrium with 



the ISM. Clarke & OeyM2002) evaluated the validity of this 



assumption by noting that in galaxies like the Milky- Way, 
the radius at which cooling becomes important is larger than 
the radius at which the super-bubble comes into pressure 
equilibrium. Of importance here is the clustering of SNe, 
which concentrates the mechanical output into small regions 
of the ISM. At high redshift, cooling is expected to be more 
efficient in the much denser ISM, although this will be offset 



by the lower gas metallicity. Mac Low & McCray ( 1988 1 
derived the cooling radius to be 



R c = 350pc 



10 38 erg/s 



4/11 



7/11 a-27/22 



(14) 



where L is the mechanical luminosity of the SNe cluster and 
£ = Z/Zq is the ISM metallicity in units of the solar value. 
Taking L = N e E a n/t e , 



R c = 0.12 kpc 



11 / t 



10 J V 1051er 87 V 4 x 10 V 

\ \ Ti x 14 / 14 

A \ 11 I TTld \ ~TT / C s 



0.05/ V0.17/ V 10km /s 

( M \ 33 (l + z\ 11 ( C 
\10 S M Q ) \~L0~ 



0.05 



(15) 



Comparing equation ( 15 1 with equation ( 13 1 for R c we find 



that the assumption of adiabatic expansion is valid in the 



low mass galaxies thought to drive reionization. We therefore 



adopt equation ( 13 1 for R c in the remainder of this paper. 



In the limit where SNe evacuated regions are smaller 
than the scale height of the disk, and the starburst lifetime 
tsF is much larger than the gas evacuation timescale i e , the 
fraction Fsn of the SNe energy may be used in feedback 
suppressing subsequent star formation. However, if the SNe 
evacuated regions break out of the disk, or tsp < t e , not 
all of the energy will be available for feedback. Based on 
the ISM porosity model of |Clarke fc Oey] ( |2002[ ) , a fraction 
fa = 2H/R e of the SNe energy goes to increasing the ISM 
porosity for disks where R c > H. In this case we find 



fa = 0.85 



Nc 
10 

M 
10 8 M G 



/ Est 



V 10 51 erg 
l + z 



A 
O05 



ma 
0.17 



10 



V10km/s 



(16) 



as long as fa < 1 and fa = 1 otherwise. Similarly, in cases 
where £sf < t e ~ 4 x 10 7 yrs, only 



f t = (t SF /t e y 



(17) 



of the overall SNe energy output is generated by the time 
the starburst concludes. The quadratic dependence on time 
arises because the number of bubbles produced grows in 
proportion to time, while the maximum size of a bubble at 



time t < t c is also proportional to time ( Oey & Clarke 1997 1. 
In cases where £sf > t e we have ft = 1. 

In the next section we fit our model to the recent data 



of Smit et al. (20121 in order to constrain tsF and /* 



3 RESULTS 

Observations of the star formation rate density function 
(SFRD) are plotted in Figure [l] for four different redshifts 
z ~ 4, 5, 6 and 7. From these observations we constrain the 
two free parameters of our star formation model, /*, ma x and 
tsF- We use the model to calculate SFRD functions for com- 
binations of these parameters and calculate the x 2 °f the 
model as 



X (/*,max, 
N, 



E 



tsF) = 

log* (SFIU,f* 



,t SF ,z) - log t> ohB {SFRi,z) 



asF R (SFRi,z) 



(18) 



Here 't> oc , B (SFRi, /*, m ax, £sf, z) is the observed star forma- 
tion rate density measured at redshift z, with uncertainty 
in dex of asFR(SFRi). In calculating likelihoods at z ~ 4 
and z ~ 5 we increased the quoted error bars by factors 
of 3 and 2 respectively in order to obtain a reduced \ 2 °f 
order unity. The SFRD function is sensitive to the value of 
Fsn, and we therefore integrate the likelihood over a range 
of values uniformly distributed between — 1 < log 10 Fsn < 



,tsp) oc 



d (l°Sio Fsn)e 



-x 2 /2 



(19) 



We note that the relation between SFR and M in equa- 
tion |l| is not perfect. As part of our comparison with ob- 
servations, and to account for scatter in this relationship, we 
convolve the predicted SFRD function equation |2| with a 
Gaussian of width 0.5 dex in SFR. An intrinsic scatter of 0.5 
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Figure 1. Comparison between the observed and modelled SFRD function (plotted as \t = In 10 X SFR X <E>). The four panels show 
results for different redshift values. In each panel, the four curves correspond to a different choice of model parameters *sp and /* imax , 
labeled by the symbols in Figure [2] The thick solid lines represent values close to the best fit. 



dex is motivated by the scatter in stellar mass at constant 
SFR found by |Gonzalez et aT"| ( |201l| . However, in addition 
we find that the value of 0.5 dex provides the best statis- 
tical fit to the observations. Our qualitative results are not 
sensitive to the choice of this scatter, however. 



3.1 Parameter constraints 

Constraints on /*, max and tsF for this model are shown in 
Figure [2] We note that these constraints are formal values 
assuming the particular model (for example, specifying the 
scale of a major merger). They are therefore indicative of the 
parameters describing the star formation model, whereas the 
overall uncertainty on starburst lifetime and efficiency may 
be larger than the ranges shown. We find that the shape 
of the SFRD function requires starburst durations of a few 
tens of Myr at z ~ 5, z ~ 6 and z ~ 7, with a few per- 
cent of the gas turned into stars per burst. For compari- 
son the left and right hand vertical grey regions represent 
times smaller than the lifetime of the most massive stars 



(ts ~ 3 x 10 6 years), and times in excess of the lifetime of 
the least massive stars that produce SNe respectively. Our 
results therefore indicate that star formation in high redshift 
galaxies is terminated on the same timescale as feedback 
from SNe can be produced. The constraints can be under- 
stood qualitatively. First, larger values of /*, max at fixed tsv 
lead to smaller values of halo mass at fixed SFR, and hence 
a larger SFRD function. Similarly, a longer starburst lifetime 
and hence duty cycle, implies that smaller halo masses are 
needed for a given value of SFR. Figure [I] shows the com- 
parison between observed and modelled SFRD functions. 
The four curves shown correspond to model parameters tsF 
and /*,max labeled by the symbols in Figure [2] We show the 
case of Fsn = 0.3, except at z ~ 4 where the lower value of 
Fsn — 0.1 produces a better fit. The thick solid lines provide 
the best fit to the observational data. The agreement with 
the data at multiple redshifts (z = 5, 6 and 7) is impressive 
since the values of /*, m ax and tsF were kept fixed. We find 
that the SFRD function should continue to increase to much 
fainter levels than currently observed. 
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Figure 2. Constraints on the model parameters /* jma x and tgF at four different redshifts (constraints are independent at each redshift). 
In each case, three contours are shown corresponding to differences in x 2 relative to the best-fitting model of A\ 2 = X 2 — Xmm = ^> 2.71 
and 6.63. Projections of these contours on to the axes provide the 68.3, 90 and 99 per cent confidence intervals on individual parameter 
values. The vertical grey regions represent time-scales longer/shorter than the lifetime of the highest/lowest mass SNc progenitor (3 X 
10 6 yr/4 x 10 7 yr). 



The parameters /*, ma x and £sf refer to single bursts, 
whereas our model includes multiple bursts at the rate of 
major mergers. We therefore calculate the total star forma- 



et al. 2010 Mufioz & Loeb 20111. Halos thought to host 



tion efficiency /*,t 



N a 



e /* (i.e. the sum of /* over 



all mergers), as well as the overall duty-cycle iV m erge£sF/tH- 
These quantities are plotted in Figure [3] based on our model 
with parameter choices corresponding to the examples in 
Figure [T] We find duty-cycles of a few percent, with higher 
duty-cycles at higher redshift reflecting the increased ratio 
between the lifetime of massive stars and the age of the Uni- 
verse. We find that ~ 5— 10% of the gas forms stars in bright 
galaxies of SFR ~ 1 — lOOM© per year, with lower fractions 
down to a percent in fainter galaxies. 



3.2 The halo masses of star- forming galaxies 

In Figure [1] we show the relation between SFR and halo 
mass M for galaxies based on our model with parameter 
choices corresponding to the examples in Figure [I] Galaxies 
observed to have SFR ~ lM0/yr reside in halos o f M ~ 
10 10 M Q , a result consistent with previous work (e.g. |Trenti 



the smaller galaxies at the hydrogen cooling limit (M ~ 
1O 75 M ) are predicted to have low star formation rates of 
SFR ~ 10- 3 M Q /yr. 



3.3 The stellar masses of star-forming galaxies 

In Figure [5] we show the relation between SFR and stel- 
lar mass M+ = md/*,totM for galaxies based on our model 
with parameter choices corresponding to the examples in 
Figure [l] SFR ~ lMo/yr characterise galaxies with M* ~ 
1O 8 M0. Our model yields a close to linear relation between 
stellar mass and star formation rate, in good agreement with 
simulations (e.g. |Jaacks et al.|2012 Salvaterra et al.|2011 l. 

In Figure [5] we also show data points representing the 
mean relation between observed star formation rates and 
stellar masses. The relation between stellar masses and 
extinction-uncorrected star formation rates at z ~ 4, 5 and 
6, was presented in Gonzalez et al. (20111. We have cor- 



rected these star formation rates for extinction using the 



methods outlined in Smit et al. (2012 L based on relations in 
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Figure 3. The values of total star formation efficiency /*,tot (i.e. the sum of /* over all mergers), and the overall duty-cycle (i.e. the 
fraction of a Hubble time during which a galaxy is starbursting) as a function of SFR. The four curves shown correspond to the SFRD 
functions shown in Figurejl] with model parameters rgp and /* jma x designated by the symbols in Figure[2] 



Meurer et al. ( 1999 1 and Bouwens et al. ( 2012 I. We find that 



our predicted SFR — M+ relation is in agreement with the 
observations. In particular, we note that while the relation 
is fairly insensitive to the value of /*, m ax, larger values of 
tsF than suggested by our modelling of the SFRD function 
would imply a star formation rate at fixed stellar mass that 
is lower than required by the observations. 

In Figure [6] we show the relations between the halo 
mass and stellar mass. The stellar to halo mass ratio is con- 
stant at high masses (M ;>1O 9 M0), but steeper towards low 
masses owing to the SNe feedback lowering the star forma- 

( [2012] ) 



tion efficiency in the model. Recently, Rocha et al 



suggested that the Milky- Way dwarf spheroidals appear to 
have had their star formation quenched at the time of infall 
into the Milky- Way, and that this time of infall was typi- 
cally ~ 7 — 10 Gyr ago. If true this implies that the Milky- 
Way dwarf spheroidals represent fossil records of the star- 
forming galaxies during reionization, and should have stellar 
masses described by our model. To make this comparison 
we therefore over-plot the observed relation for Milky- Way 
dwarf spheroidals ( Boylan-Kolchin et al. [[2012 1 . To convert 
to halo mass we adopt the values of virial velocity at time 
of collapse from Boylan-Kolchin et al. (20121, and calcu- 



late the corresponding halo mass at each redshift z. We find 



good agreement between the observed relation (although the 
scatter is large) , indicating that if the dwarf spheroidals are 
old galaxies formed at around the end of reionization, then 
they would have the stellar to halo mass ratio predicted by 
our SNe feedback limited model. Combined with the cor- 
rect prediction of the stellar mass for star-forming galaxies 
during reionization, this implies that our model correctly 
describes the stellar mass to halo mass ratio in the range 
10 5 M,3 <M* <10 10 M B . 



4 THE ESCAPE FRACTION OF IONIZING 
PHOTONS 

The contribution of galaxies to reionization is dependent on 
the star formation rate and initial stellar mass function, but 
is also limited by the fraction of ionizing photons that escape 
their host galaxies. If the escape fraction is small, then star 
formation had to be very efficient at high redshift in order 
to reionize the Universe. The escape fraction is therefore a 
critical parameter in studies linking high redshift galaxy for- 
mation to reionization. In this section we discuss the impli- 
cations of our model for the ionizing photon escape fraction. 
Attempts to determine the ionizing photon escape frac- 
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Figure 4. Relation between SFR and halo mass M at four different redshifts based on our model with parameter choices corresponding 
to the examples in Figure [l] Also shown is the critical star formation rate required to achieve porosity equal to unity (5FR cr ; t ) as a 
function of halo mass. Regions of the M — SFR plane where the porosity < 1, and hence ionizing photons do not escape, are shaded 
grey. 



tion have been dominated by direct observations of rela- 
tively low redshift galaxies, and by numerical simulation. 
Observational estimates of the escape fraction at the Lyman- 
limit are currently uncertain, with no confident detections 
at < 2 < 1, and only some detections at z > 3. At red- 
shifts z ~ 1 — 3, observations have suggested a broad range 
of values for / csc , from a few percent to ;>20% (ISteidel et al.l 
2001) JFernandez-Soto et al.|2003||Shapley et al.|2006||Siana| 



et al.|2007 e.g.). Inoue et al. (20061 have examined the evo- 
lution of the escape fraction in the redshift range z — — 6 
using both direct observations of the escape fraction and 
values that they derive from measurements of the ionizing 
background. They find that the escape fraction evolves from 
/csc ~ 1 — 10%, increasing towards high redshift. 



Theoretically, Razoumov & Sommer-Larsen ( 2006 1 used 



galaxy formation simulations incorporating high-resolution 
3-D radiative transfer to show that the escape fraction 
evolves from / eS c ~ 1-2% at z = 2.39 to f esc ~ 6-10% at 



z = 3.6. In agreement with Fujita et al. ( 2003 1, Razoumov & 



Sommer-Larsen ( 2006 1 (see also Yajima et al 



2009) find that 



increased supernova feedback at higher redshift expels gas 
from the vicinity of starbursting regions, creating tunnels in 
the galaxy through which ionizing photons can escape into 
the IGM. Numerical simulations ( Gnedin et al.|[2007 1 have 
predicted a value for / eS c between 1 and 3%, for halos of mass 



M £5 x 1O 1O M , over the redshift range 3 < z < 9. This 
very low efficiency of reionization would have profound im- 
plications for the reionization history. In addition to a small 
escape fraction in massive galaxies, Gnedin et al. (20071 
further predict that halos with M <J5 x lO iu M0 have an 
escape fraction that is negligibly small. However, more re- 



cently |Wise fc Cen ( 2009 1 have used a large suite of simu- 



lations to show that the time averaged escape fraction for 
dwarf galaxies is expected to be large (> 25%). Since dwarf 
galaxies are thought to dominate the ionizing flux, resolution 
of this issue is of primary importance for studies of reion- 
ization. Overall, the escape fraction is predicted to span a 
very broad range, 0.01 jj/cac <;1- This broad range may be 
explained by inhomogeneities in the hydrogen distribution 
within galaxies ( |Dove et al.|2000||Fernandez fc Shull|2011[ ), 
or by variations in viewing angle ( |Wood fc Loeb||2000 l. 



Within our formalism we follow Clarke & Oey ( 2002 1 
who proposed that galaxies with a sufficient star formation 
rate to generate a porosity greater than unity had an escape 
fraction for ionizing photons that is of order unity, whereas 
galaxies with insufficient start formation have a negligible 
escape fraction. The critical star formation rate required to 
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Figure 5. Relation between SFR and stellar mass M* at four different redshifts based on our model with parameter choices corresponding 
to the examples in Figure]!] Also shown is the critical star formation rate required to achieve porosity of unity (SFR cr [ t ) as a function of 
halo mass. Regions of the — SFR plane where the porosity < 1, and hence ionizing photons do not escape are shaded grey. The points 
show the observed relations between stellar masses and star formation rates. We have corrected the star formation rates in Gonzalez 
for extinction using the methods outlined in |Smit et al.|(2012} , based on relations in |Meurer et all] 41999} and |Bouwens| 



et al. 


42011} 


et al. 


42012} 



achieve porosity of unity is ( Clarke fe Qey|2 002 ) 



SFR C 



O.lSMoyr- 1 (/ t / d 



m d M \ / 



10 10 M o y Vl0km/s 



(20) 

In Figures [1] and [5] we plot the curves corresponding 
to SFR CI it as a function of halo mass and stellar mass, re- 
spectively. We shade the areas below the line for clarity to 
indicate the regions in the M — SFR and M* — SFR planes 
where the porosity < 1 and hence ionizing photons do not 
escape. The model predicts that galaxies with halo masses 
in excess of M ~ 10 8 M Q (corresponding to M t > 1O 5 M ) 
have a sufficiently large SFR to potentially contribute to 
reionization. 



In the model of Clarke & Oey ( 2002 1 galaxies with 
SFR > SFR cl it attain porosity of unity after a time 
t Q = t c (SFR/SFR clit )~ 1/2 . The ionizing photon rate from 
a starburst drops rapidly once the most massive stars fade 
away on a timescale of t B ~ 3 x 10 6 years ( |Barkana fc Loeb 



2001 ). Since ionizing photons are produced for a time tsF-Ns 



but can only escape the galaxy after a time tQ, the time dur- 
ing which the galaxy with SFR > SFR CI it emits photons 



into the IGM is 



tion — tSF + t B — tQ 



: - [ tA ( SFR x 



ts J \SFRcrit 



(21) 



For small values of tss < t e , Q may never reach unity for 
some values of SFR, even if SFR > SFR cr i t - To estimate 
tion we must therefore integrate over the distribution for 
SFR at fixed halo mass, which was assumed to have a scat- 
ter of 0.5dex. Thus we calculate 



(tion) = / dAsFR ti on (SFR x 1CT 



2x0.5^ _ (22) 



We note that since (tion) < tsF + t s , we expect a wide range 
of observed escape fractions for star forming galaxies de- 
pending on whether or not the starburst is being observed 



before or after tQ . The galactic porosity model of ( Clarke & 
Oey||2002 1 postulated that the escape fraction is negligible 
at t < tQ and of order unity at t > tQ. Our model provides 
the probability Pi on that a galaxy with a particular SFR 
should be observed with a non-zero escape fraction for ion- 
izing radiation, calculated as the fraction of time that the 
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galaxy is observed to be star-forming and for which t>tq, 

Pion = ^ 
tSF + ts 

Curves of P on are plotted as a function of SFR in Figure [7] 
The left hand panel shows the best fit case of tsF = 10 7 yr, 
and the right hand panel tsF = 2 x 10 7 yr. Four redshifts are 
shown. Our model predicts that large escape fractions should 
be rare for both low SFR and high SFR galaxies. However, 
we expect approximately a half of star-forming galaxies with 
SFR ~ 0.1 — IMo/yr to have a significant escape fraction. 
Thus, our model provides a natural explanation for the wide 
range of conclusions regarding observations of the escape 
fraction of ionizing photons from star-forming galaxies. 

Next we investigate the implications of this model for 
the reionization history. 



5 IMPLICATIONS FOR THE REIONIZATION 
HISTORY 

We continue with a semi-analytic calculation of the reioniza- 
tion history of the IGM based on the star formation model 
presented in § [2] The basis for our model of reionization is 
the excess ionization rate over the recombination rate for 
hydrogen in an inhomogeneous IGM. [Miralda-Escude ~et al.| 



( 2000 ) presented a model that allows the calculation of an ef- 
fective recombination rate in an inhomogeneous universe by 
assuming a maximum overdensity (A c ) penetrated by ion- 
izing photons within HII regions. The model assumes that 
reionization progresses rapidly through islands of lower den- 
sity prior to the overlap of individual cosmological ionized 
regions. Following the overlap epoch, the remaining regions 
of high density are gradually ionized. It is therefore hypoth- 
esized that at any time, regions with gas below some critical 
over density Ai = pi/{p) are highly ionized while regions of 
higher density are not. The fraction of mass in regions with 
overdensity below Ai, is found from the integral 



i*M ( Ai) 



dAP v (A)A, 



(24) 



where Py(A) is the volume weighted probability distribution 
for A. Miralda-Escude et al. ( |2000[ ) quote a fitting function 
that provides a good fit to the volume weighted probability 
distribution for the baryon density in cosmological hydro- 
dynamical simulations. In what follows, we draw primarily 
from the prescription of |Miralda-Escude et al.| ( |2000[ ) and re- 
fer the reader to the original paper for a detailed discussion 



of its motivations and assumptions. Wyithe & Loeb (2003a I 



employed this prescription within a semi-analytic model of 
reionization. This model was extended by |Srbinovsky fc| 
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Wyithe (20071 and by Wyithe et al. (20081. We refer the 



reader to those papers for a full description. 

The quantity Qi is defined to be the volume filling factor 
within which all matter at densities below Ai has been ion- 
ized. The reionization history is quantified by the evolution 
of Qi that follows the rate equation 



dQj 

dz 



1 



dn~. 



n F M {Ai) dz 

a B (l + z) J?(Aj)n — H — - 

dz dz 



f?5) 



where q?b is the case B recombination coefficient, no is the 
mean comoving density of hydrogen in the IGM, and R(Ai) 
is the effective clumping factor of the IGM. The evolution 
is driven by the rate of emission of ionizing photons per co- 
moving volume dn~ t /dz. Within this formalism, the epoch 
of overlap is precisely defined as the time when Qi reaches 
unity. Prior to the overlap epoch we must solve for both Qi 
and Fm (or equivalently Ai). The relative growth of these 
depends on the luminosity function and spatial distribution 
of the sources. In this regime we assume Ai to be constant 
with redshift before the overlap epoch and compute results 
for models with values of Ai = A c = 10. Different values 
of A c are not found to quantitatively effect our results for 
values in the range 5 < A c < 20 ( |Wyithe et al.|2008[ ). 

Following overlap we may describe the post-overlap evo- 
lution of the IGM by computing the evolution of the ionized 
mass fraction according to the equation 



dF M ( Aj) 
dz 



1 dn~f 
no dz 



-a B (l + z) 3 R(Ai)n ^. 



(26) 



This follows directly from equation (25 1 with Qi — 1. In this 
post overlap regime the value of Ai is the dependent variable 
describing the ionization state of the IGM (whereas prior to 
overlap Ai = A c ). Equation (261 is integrated to obtain Fm 
(or equivalently Ai) as a function of redshift. 

The emission rate of ionizing photons per co-moving 
volume that is required to compute the reionization history 
can be written 



dn 7 
dz 



N. 



(/esc/?*) dt 

m D dz ' 



(27) 



where iV 7 is the number of ionizing photons produced per 



baryon incorporated into stars, and (/ esc p*) is the product 
of the escape fraction and the star formation rate density 
averaged over halo mass. As described in the introduction, 
only a fraction / esc of the ionizing photons produced by stars 
enter the IGM. We define (p* lCSC ) using our formalism as 



r 3o 
J Mi 



dM /esc SFR x 



in 



dN 2 



dtdAM 



Mi, AM 



dn 
dM 



where (ti 0n ) is the time during which the galaxy emits ion- 
izing photons into the IGM (equation 22 1 . 

In a cold neutral IGM beyond the redshift of reioniza- 
tion, the collapsed fraction should be computed for halos 
of sufficient mass to initiate star formation. The minimum 



virial temperature is set by the temperature T m 



10 4 K, 



above which atomic hydrogen cooling promotes star forma- 
tion. Following the reionization of a region, the Jeans mass in 
the heated IGM limits accretion to halos above Ti on ~ 10 5 K 



( Efstathiou||1992| |Thoul fc Weinberg||l996| |Dijkstra et aL 
|2004[ ). Including these separate components from ionized 
and neutral regions of IGM, we get 

(/cscp*) = </e S c/5*)[>M mln ](l-<2i) + (/esc/9*) [>M ian ]Qi, (28) 

Our model assumes a spectral energy 
Population-II stars with a metallicity Z - 



distribution of 
= O.O5i?0 and a 
Scalo ( 1998 1 initial mass function, for which the resulting 



number of hydrogen ionizing photons per baryon incorpo- 
rated into stars is iV 7 ~ 4000 (Bar kana fc Loeb||2001 1. 

In order to estimate the ionizing background following 
the end of reionization, we compute a reionization history 
given a particular value of A c , combined with assumed val- 
ues for /* and / csc . Given this history, we then compute the 
evolution of the background radiation field due to the same 
sources. After the overlap epoch, ionizing photons will expe- 
rience attenuation due to residual over dense pockets of HI 



gas. We use the prescription of Miralda-Escude et al. ( 2000 1 



to estimate the ionizing photon mean- free-path. Following 
Oh & Furlanetto (20051, who note that the the constant of 



proportionality relating mean-free-path to the volume- filling 
fraction should be reduced by approximately a factor of 2, 
we adopt A = (30kms -1 )/ff(l — Q)~ 2 ^ 3 , and subsequently 
derive the attenuation of ionizing photons. We then compute 
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Figure 8. Models for the reionization of the IGM and the subsequent post-overlap evolution of the ionizing radiation field. Four cases 
are shown, with line styles corresponding to the cases over plotted on the constraints on / tjmax and t^p at z ~ 6 that are repeated in the 
top left panel for reference. A value of A c = 10 was adopted for the critical overdensity prior to the overlap epoch. Upper Right Panel: 
the evolution of volume averaged (dark lines) and mass averaged (grey lines) neutral fraction. Lower L eft Panel: Ionization rate as a 
function of redshift (in units of 10 _12 s _1 ). The observational points are from |Bolton fc Haehnelt| j2007| and|Wyithe fc Bolton| (201l| ). 
Lower Right Panel: The mean-free-path for ionizing photons. The data point is based on|Storrie-Lombardi et al.| l|1994[|. 



the flux at the Lyman-limit in the IGM due to sources imme- 
diate to each epoch, in addition to redshifted contributions 
from earlier epochs. 

We show the results of this modelling in Figure [8] for 
four models, including the best-fit case presented in Fig- 
ure [T] For reference the constraints on /*. max and tsF at 
z ~ 6 are repeated in the top left panel. The upper right 
panel shows the evolution of volume averaged (dark lines) 
and mass averaged (grey lines) neutral hydrogen fraction. 
To illustrate that the model for star formation is consistent 
with the reionization history, we also show the observables of 
ionization rate as a function of redshift (lower left), and the 
mean-free-path for ionizing photons (lower right). We find 
for a choice of / csc = 0.8, that a model near the best fit for 
the SFRD function (tsF = 2 x 10 7 yr) results in a reioniza- 
tion history that is consistent with low redshift constraints 
(thick lines). The reionization history shown is consistent 
with recent constraints from the patchy kinetic Sunyaev- 
Zel'dovich effect ( Zahn et al.|2 012 ). The model produces an 
optical depth to electron scattering of r os = 0.065, which is 
lower than the observed value of r cs = 0.088 ± 0.015 (Ko- 



matsu et al.|2011 1. However, our model does not include the 
possibilities of population-Ill stars or an increased escape 
fractions at very high redshift ( jAlvarez et al.|[2012[ ). Mod- 
els with larger /*, max reionize the universe too early, while 
those with smaller /*, max or larger £sf reionize the universe 
too late. For the model with a smaller value of tsF = 10 7 yr, 
we find that reionization occurs very late (at z ~ 5), and 
greatly underestimates the ionizing background. This is be- 
cause porosity of unity is rarely achieved by the time the 
starburst is finished in this model, so that very few ionizing 
photons escape. This indicates that the reionization of the 
IGM may limit the starburst time-scale to be 4sf J>10 7 yr. 

In the left panel of Figure [9] we show the SFR 
weighted average of the escape fraction (i.e. {f aS c)sFR. = 
{fescp*)/{p*}) as a function of redshift for the four mod- 
els shown in Figure [8] Since the value of escape fraction in 
highly star forming galaxies is / esc = 0.8 in this model, we 
find that most ionizing photons produced in galaxies are be- 
ing lost either because the SFR was insufficient to attain 
porosity equal to unity, or because the UV luminous stars 
had died before porosity of unity was achieved. Interestingly, 
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Figure 9. Contribution of galaxies to rcionization. Left Panel: The evolution of SFR weighted average escape fraction for the four 
models shown in Figure [8] Right Panel: The evolution of volume averaged neutral fraction. Three cases are shown, including the best 
fit model from Figure [8] (thick black line) for comparison with models in which photons are allowed to escape from all galaxies, with 
/ csc = 0.8 (i.e. SFR clit = 0). In these cases we show examples where radiative feedback is included (fion — 30km/s and Vi n — 

50km/s), 

and a case in which it is not (i.e. «; on = 0). 



we find that the SFR averaged escape fraction increases to- 
wards high redshift, being twice as large at z ~ 10 as at 
z ~ 4. This trend is consistent both with the requirements 
of the Lya forest at z ~ 4 (Bolton & Haehnclt 2007J, and 
recent empirical estimates for the escape fraction of ioniz- 
ing photons from high redshift galaxies during reionization 
( Finkelstein et al.|2012 i. 



THE CONTRIBUTION OF LOW MASS 
GALAXIES TO REIONIZATION 



A range of models (e.g. Haiman fc Holder||2003 Iliev et al 



2007) have predicted that radiative feedback, which can sup- 



press star formation in very low mass galaxies, will have a 
significant effect on the reionization history. The effect is 
often referred to as self regulation of reionization, because 
it could lead to the reionization process being extended. 
In these models reionization is predicted to start at high 
redshift due to very low mass galaxies with virial temper- 
atures of ~ 10 4 K. However, once regions become ionized 
subsequent star formation in these very low mass galaxies is 
quenched, so that reionization must be completed by galax- 
ies with a virial temperature ;>10 K which form later. The 
reason why radiative feedback could have a large effect on 
the reionization history is well understood. At high redshift 
the collapsed fraction of dark matter halos is dominated by 
low mass galaxies with a virial temperatures smaller than 
~ 10 5 K. Despite the focus of numerical and semi-analytic 
studies of high redshift galaxy formation on the role of feed- 



back on star formation in low mass galaxies (e.g. Raicevic 



et al. 2011 Finlator et al. 20111, models of the structure 



of reionization have generally assumed the luminosity to be 
proportional to halo mass, and that the escape fraction of 



ionizing photons is independent of mass (e.g. Iliev et al. 2007 
McQuinn et al.|2007 l. As a result, the ionizing photon bud- 



get for reionization in these constant mass-to-light models 
is also dominated by low mass galaxies prior to reioniza- 
tion. However, in this paper we have found that SNe feed- 
back is required to reproduce the shape of the high redshift 



SFRD function, in agreement with the previous numerical 
work. Thus, the contribution to reionization from low mass 
galaxies may be smaller than previously thought ( | Raicevic] 
et al. 112011 1. In this section, we investigate the contribution 



to reionization from galaxies of various masses. 

The effect of low mass galaxies on the reionization his- 
tory is illustrated in the right panel of Figure [9] which shows 
the evolution of volume averaged neutral fraction. The best 
fit model from Figure [8] is shown (thick black line) for com- 
parison with models in which the critical SFR criteria for 
ionizing photons to escape is not applied so that photons 
are allowed to escape from all galaxies with / osc = 0.8 (grey 
lines). For the latter, three cases are shown including where 
radiative feedback is not included (i.e. Vi on = 0) and where 
"ion = 30km/s and fJi on = 50km/s. We see that the appli- 
cation of the critical SFR criteria has a significant effect on 
the reionization history, leading to an earlier reionization by 
Az ~ 1.5 at a fixed / csc . However, suppression of galaxies 
smaller than Vi on = 30km/s or 50km/s results in a smaller 
change to the reionization history relative to a model in 
which radiative feedback does not operate (Az <^0.2), a find- 
ing that is contrary to the standard picture for self-regulated 
reionization (e.g. Iliev et al.|2007[ ) . 

To understand the reasons for this we show the cumu- 
lative contribution to the number of ionizing photons per 
hydrogen in the Universe per Hubble time as a function of 



halo virial velocity (left panel of Figure 10 1. In calculat- 
ing these curves we have not applied the SFR cr i t criteria 
(equation 20 1, and assumed that ionizing photons escape all 
galaxies with / csc . Contributions at three different redshifts 
are shown. In the right panel we repeat these results, but 
plot the contributions as a function of star formation rate 
rather than virial velocity. We find that photons from galax- 
ies below Uvir ~ 30km/s (with SFR < 0.01M Q /yr), corre- 
sponding to those affected by radiative feedback, represent 
only ~ 10% of the potential photon budget at z ~ 6 and 
only ~ 30% at z ~ 10. Thus, these low mass galaxies make 
only a small contribution to reionization, implying that ra- 
diative feedback should not be important in regulating the 



© 0000 RAS, MNRAS 000, 000-000 



14 Wyithe & Loeb 




Figure 10. Contribution of low mass galaxies to reionization. Left Panel: The cumulative contribution to the number of ionizing 
photons per hydrogen in the Universe as a function of halo virial velocity. To investigate the possible contribution to reionization from 
small galaxies we have not applied the SFR crl t criteria, and instead assumed that ionizing photons escape all galaxies with / osc = 0.8. 
Contributions at three different redshifts are shown, and expressed in units of photons per hydrogen in the IGM times f^ (in order to 
make the curves independent of / CS c)- Right Panel: Results repeated as a function of the SFR per galaxy. 



reionization process. This reduced contribution arises be- 
cause SNe have lowered the star formation efficiency to a 
level where galaxies with u v h- ~ lOkm/s contribute a very 
small fraction of the total star formation rate, despite the 
corresponding halos representing the dominant component 
of collapsed fraction of dark matter. 

Thus, we conclude that SNe feedback was much more 
important than radiative feedback in shaping the reioniza- 
tion history, a result that is consistent with both |Raicevi6 



Kim 



et al. (2011 1 and the recent semi-analytic modelling of 
et al.| ( |20r2| ) using GALFORM within the Millennium-II 
simulation. We note that an important caveat to these 
results is the possibility of a top-heavy mass function of 
Population-Ill stars in small galaxies. Since in this case 
the ionizing efficiency is much larger than for Population-II 
stars, Population-Ill stars in early low-mass galaxies could 
make a more important contribution to reionization, par- 



tially reionizing the Universe at earlier times (e.g. Wyithe 
fc Loeb|[2003a| |Cen||2003| >. Overall we find that half of the 



ionizing photons are produced by galaxies with v v i r J>60km/s 
corresponding to SFR J>lM0/yr. These are observed galax- 
ies, indicating that observations summarised in |Smit et al.| 
( 2012 1 correspond to about half of the ionizing photons pro- 



duced by galaxies. This result is consistent with the recent 



findings of Finkelstein et al. (20121, who suggested that the 



observed population of z ~ 6 galaxies is sufficient to pro- 
vide most of the ionizing photons required for reionization 
assuming escape factions below 50%. 



7 CONCLUSION 

We have described a simple model for the star formation 
rate density function at high redshifts based on the extended 
Press- Schechter formalism. This model assumes that a star- 
burst is triggered by each major merger, lasting for a time 
tsF and converting at most /* jmax of galactic gas into stars. 
We include a simple physical prescription for SNe feedback 



that suppresses star formation in low mass galaxies, and re- 
sults in a minimum galaxy mass from which ionizing photons 
can escape. This model for star formation has only two free 
parameters, tsF and /V.max, but accurately describes recent 
measurements of the amplitude and shape of the star for- 
mation rate density function between redshifts 4 and 7. 

Comparison of our modelling with observational data 
implies that individual starbursts were terminated after 
tsF ~2x 10 7 years. This termination time lies between the 
main-sequence lifetimes of the lowest and highest mass SNe 
progenitors, indicating that starbursts were quenched once 
SNe feedback had time to develop. From the sum of all ma- 
jor merger events, high redshift galaxies convert ~ 5 — 10% 
of their mass into stars for large galaxies with star formation 
rates above ~ lM0/yr. However, the overall star formation 
efficiency of lower luminosity galaxies is only a few percent. 
In our model, high redshift galaxies have a high duty-cycle 
for star formation, undergoing starbursts ~ 10% of the time 
atz~ 10. 

We calculate the relation between stellar mass and star 
formation rate, finding it to be approximately linear, in 
agreement with previous theoretical work and observations. 
The predicted ratio of star formation rate to stellar mass 
based on our fit to the star formation rate density function 
agrees with the observed values for star-forming galaxies at 
z ~ 4-6 where stellar mass has been measured. Moreover, if 
the Milky- Way dwarf spheroidals represent fossil records of 



the star-forming galaxies during reionization (Rocha et al 



2012 1, they should also have stellar masses described by our 



based on the galactic porosity model of Clarke & Oey ( 2002 1 



model and we find good agreement with the predicted rela- 
tion between stellar and halo masses (although the scatter 
is large). This implies that our model correctly describes the 
stellar mass to halo mass ratio from 10 5 M Q <CM* <TO 1O M . 

We used our model to discuss the escape fraction of 
ionizing photons from high redshift star-forming galaxies. 
|Clarke Sz Oey| ( |2002[ ) proposed a simple model where galax- 
ies with a sufficient star formation rate to generate a poros- 
ity greater than unity had an escape fraction for ionizing 
photons that is of order unity, whereas galaxies with insuffi- 
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cient start formation have a negligible escape fraction. Since 
a fraction of the galaxies starburst lifetime is required to 
build the porosity to unity, every star-forming galaxy has 
a non-zero probability of being observed with a zero escape 
fraction for ionizing radiation. Our model predicts that large 
escape fractions should be rare for both low SFR and high 
SFR galaxies. However, we expect approximately a half of 
star-forming galaxies with SFR ~ 0.1 — lMo/yr to have 
significant escape fraction. Thus our model provides a nat- 
ural explanation for the wide range of conclusions regard- 
ing observations of the escape fraction of ionizing photons 
from star-forming galaxies ( jSteidel eF al. 2001 
Soto et al.|2003||Shapley et al.|2006||Siana et at 



Fernandez- 1 



2007 



e-g-J 



and predicts that the escape fraction during reionization at 
2 ~ 10 was twice as large as at z ~ 4. 

We also used a semi-analytic model for the reionization 
process based on our SNe regulated star formation history. 
Even after allowing for the fact that SNe feedback enables 
escape of ionizing radiation only after the galactic porosity 
reaches unity, we find that our model is able to reionize the 
Universe within current observational constraints. We find 
that low mass galaxies were minor contributors to reion- 
ization, owing to the suppression of star formation by SNe 
feedback. We further find that this SNe feedback lowers the 
efficiency of star formation in low mass galaxies to such an 
extent that photo-ionization feedback on low mass galaxy 
formation does not significantly effect the reionization his- 
tory. This is because the galaxies that would have been sub- 
ject to radiative feedback are only very minor contributors 
to the potential ionizing photon budget once SNe feedback 
is taken into account. F] Finally, we find that approximately 
half of the ionizing photons needed to complete reioniza- 
tion have already been observed in star-forming galaxies at 
z = 6-10. 
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